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Abstract 

We propose and test a renormalization procedure which acts in 
Hilbert space. We test its efficiency on strongly correlated quantum 
spin systems by working out and analyzing the low-energy spectral 
properties of frustrated quantum spin systems in different parts of the 
phase diagram and in the neighbourhood of quantum critical points. 
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Introduction. 



Microscopic many-body quantum systems are often subject to strong in- 
teractions which act between their constituents. Perturbation treatments 
make sense if it is possible to introduce a mean-field which is able to absorb 
the main part of the interaction between the constituents. But often such an 
approach does not lead to sensible results, in particular in the case of realis- 
tic quantum spin systems. Non-perturbative procedures are needed, see f.i. 

Spectral properties of quantum systems are obtained by means of the di- 
agonalization of a many-body Hamiltonian in Hilbert space which is spanned 
by a complete, in general infinite or at least a very large basis of states al- 
though the information of interest is restricted to the knowledge of a few low 
energy states. In order to reduce the problem to these states we proposed a 
new approach which acts as a size reduction of Hilbert space. The procedure 
consists of an algorithm which eliminates states by means of a step by step 
projection procedure to subspaces of the original space [2]. It rehes on the 
general renormalization concept [U [El [T51 [TB] . It can be applied to all types 
of microscopic quantum systems (molecular, atomic, nuclear, solid state,...) 
in contradistinction with more specific procedures such as f.i. the Density 
Matrix Renormalization Group {DMRG) [111 [I2l [ISl [7] which are more 
specifically applied to systems which live on a lattice. The reduction concept 
shows some connection with the method which has been proposed and ap- 
plied in problems related to the treatment of phonons [SI E] by means of an 
efficient method which connects a density matrix approach with optimized 
space dimension reduction. 

In the present work we implement this algorithm as a preliminary test 
of the practical efficiency and accuracy of the method when applied to large 
strongly interacting systems. We want to test its capacity to deliver precise 
information about the low energy eigenstates of the many-body systems it is 
aimed to describe. Here we use quantum spin systems as test probes. 

The reduction procedure. 

We consider a system described by a Hamiltonian depending on a unique 
coupling strength g which can be written as a sum of two terms 

H = Ho + gH, (1) 

The Hilbert space T-C^^^ of dimension is spanned by a priori arbitrary 
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complete set of basis states i = 1, - ■ ■ , N} which may be f.i. the eigen- 

states of Hq. Eigenvectors {"^f^^) decompose on the basis as 



N 

^r)=E«r(^^"^)i^^) (2) 



1=1 



where the amphtudes {0'^ji\g^^^)} depend on the value g'^^'' oi g in the chosen 
Hilbert space H'^^l 

The space can be decomposed into two subspaces by means of the pro- 
jection operators P and Q [17J, 



The projected eigenvector P|\E'^-^'') obeys then the effective Schrodinger equa- 
tion 



ff./,(Ary Iff) = Af P|*f' ) (4) 

where Heff{\^^^) is a new Hamiltonian which operates in the subspace 
Pl-L^^\ It depends on the eigenvalue A^^^ which is also the eigenenergy 
corresponding to I^E'j^'*) in the initial space T-i^^\ The coupling g^^^ which 
characterizes the Hamiltonian H'^^^ in H^^-* is now constrained to change 
into g^^~^^ in such a way that the eigenvalue in the new space 7Y(^~^) is the 
same as the one in the original space 

Af = Af ) (5) 

The determination of g^^^^^ by means of the constraint expressed by Eq. ([5]) 
is the central point of the procedure. It corresponds to a renormalization 
procedure which is implemented by means a non-linear relation between g^'^^ 
and(7(^-i) [2^,^. 

In the sequel P|\E'[^'') is chosen to be projection of the ground state eigen- 
vector |\E'[^'*) and A^^'' = A|^^~^'' = Ai the corresponding ground state energy. 
Following this first operation the reduction procedure is iterated in a step by 
step decrease of the dimensions of Hilbert space, N^-^N — l^N — 2\-^... 
leading at each step to a new coupling strength g^^^''^ which is fixed as 
the solution of an algebraic equation and is shown to obey a flow equation 
in the limit of a continuum description [21 E] . 
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The implementation of the reduction algorithm 
The procedure goes along the following steps: 

1— Consider a quantum system described by an Hamiltonian H^^^ and 
compute its matrix elements in a definite basis of states -i = 1, . . . , A^}. 

2— Use the Lanczos technique to determine A^^"* and the amplitudes 
a^^\g^'^^) corresponding to |\E'[^''((7^^-')) [19]. The diagonal matrix elements 
{ei = (<l>j I if^ I are arranged in decreasing order of values of the \aii\g''^^)\ 

3— Fix g^^~^^ as described above and in [2]. Construct H^^~^^ = Hq + 
g(^~^)Hi by elimination of the matrix elements of H^^^ involving the state 

|$7V). 

4— Repeat the procedures 2 and 3 by fixing at each step k X^-^~''^ = 
X[^^ = Ai. The iterations are stopped at = Nmin which may correspond 
to the limit of space dimensions for which the spectrum gets unstable. 

This is due to the fact that which is the eigenvector in the 

space H^^-''-^) and the projected state P\'^[^~''^) of l^j^"''^) into U^^-^-^) 
may not coincide exactly. As a consequence it may not be possible to keep 
^(Af-fc-i) j.jgQj.Qygiy equal to A^^"'^'* = Ai. In practice the degree of accuracy 

depends on the relative size of the eliminated amplitudes 

This point will be tested by means of numerical estimations and further 

discussed below. 

Quantum Phase transitions and fixed points 

Strongly correlated systems often possess rich phase diagrams and critical 
transition points EH EZ] • We show now how our procedure refiects the 
presence of these transitions. 

The eigenvalues \k{g) of H{g) = Hq + gHi are analytic functions of g 
which may show algebraic singularities [201 HH [22] at so called exceptional 
points g = ge- Exceptional points are first order branch points in the complex 
g - plane which appear when two (or more) eigenvalues get degenerate. This 
can happen if g takes values such that = ei where = In a 

finite Hilbert space the degeneracy appears as an avoided crossing for real 
g. If an energy level belonging to the P7i subspace defined above crosses 
an energy level e/ lying in the complementary QTi subspace the perturbation 
development constructed from Heff{E) diverges [22] ■ Physical states can get 
degenerate in energy for real values of g^- 
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Exceptional points are defined as the solutions of 



f{\{g,)) = det[H{g,) - X{g,)I] = (6) 

and 

^3^Ia=a(,.) = (7) 

They are fixed points of the coupling strength g which stays constant 
during the space reduction process [2l|3]. Indeed, if {Xi{g)} are the set of 
eigenvalues the secular equation can be written as 

N 

n (A - A,) = (8) 

i=l 

Consider X = Xp which satisfies Eq. ([6]). Eq. ([7]) can only be satisfied if 
there exists another eigenvalue Xg = Xp, hence if a degeneracy appears in the 
spectrum. This is the case at an exceptional point 



If the eigenvalue ''\k = 0,1,... which is either constant or con- 
strained to take a fixed value Xj gets degenerate with some other eigenvalue 
Ai^ ''\g = ge) in the space reduction process this eigenvalue must satisfy 

Ar-'='fe) = Af -"(«:) (9) 

at any step k and / of the projection procedure. In the continuum limit for 
large values of N, {N, N — 1) ^ {x,x — dx), 

dXj dXi{qJx),x) , , 

— ^ = = ^ ,^ ' 10 
dx dx 

where 

Xi{ge{x),x) = {^i{ge{x),x)\H{g,{x))\^i{ge{x),x)) 
Consequently 

dXi dXi dge ^ 
dx dge dx 

Since Aj((?e(x),x) stays constant with x in the space dimension interval 
(x — dx, x) the first term in Eq. (fTTj) is equal to zero. Hence in general 

dXi , dOf. , . 

-^^0 and -f = 12 
oge dx 
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I2 22 32 42 52 ••• h 4 i 

Figure 1: The frustrated spin ladder. The couphng strengths are indicated 
as given in the text. 



which shows that the exceptional point g = Qe corresponds to a fixed 
point in the renormalization process, and characterizes the existence of a 
quantum phase transition. 

This is general and verified for any Hamiltonian for which state degen- 
eracies occurs. Numerical applications corresponding to a linear dependence 
on g, H{g) = Hq + gHi follow. 

Applications to a frustrated two-leg quantum spin ladder. 

Quantum spin systems are strongly interacting systems whose properties 
cannot be studied by means of perturbation expansions. Hence they are 
particularly well adapted to non-perturbative treatments such as renormal- 
ization procedures. In the following we present a test of the present method 
on such a system. 

Consider a spin-1/2 ladder as shown in Figure [H Its Hamiltonian H{g) = 
Hq + gHi corresponds to Hq = 0, g = g^^^ = Jt and 

L 

Hi = ^ Si^Si^ + 'jti 5^ (si^Sj^ + Si^Sj^) + 7c isiiSj2 + Si^Sn) (13) 

i=l <ij> <ij> 

where = Ji/ Jt, Ic = Jc/ Jt which are kept constant and ^/(^) is the renor- 
malizable coupling strength. The number of sites along a leg is L. Indices 
1 or 2 label the vector spin 1/2 operators Sj^. acting on the sites i of leg /c, 
< ij > labels nearest neighbours on a leg, (ij) diagonal interactions between 
sites located on different legs. The coupling strengths Jt, Ji, Jc are positive. 
One can show that the renormalization process may be indifferently realized 
on Ji or Jc and leads to the same result [3]. 
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The basis states {|$a;), k = 1,...,A^} are generated as product states 
|l/2 m,) 

2L 

= |l/2 mi;.. .,1/2 rm; ...,1/2 m2L,^mi = Mtot) 

1=1 

with {rrii = +1/2, —1/2}. The total projection is a good quantum number 
fixed to Mtot = in the following. 

Test observables. 

In order to estimate the accuracy of the reduction procedure we introduce 
the quantity 

/ (N) _ 

p{i) = \S2i — X 100 with i = l,...,3 (14) 

where e^"'"' = /2L with n = {N — k) corresponds to the energy per site at 
the ith physical state starting from the ground state at the kth iteration in 
Hilbert space. This quantity provides a percentage of loss of accuracy of the 
eigenenergies in the different reduced spaces. 

A global characterization of the ground state wavefunction can also be 
given by the entropy per site in a space of dimension n 

1 " 

s('^) = -^5^P./rzP, wzth P,= |($,|vpS'^))|2 = |a;:)|2 (15) 

i=l 

which works as a measure of the distribution of the amplitudes in the 

physical ground state [23] . 

Results and discussion. 

We apply the reduction algorithm to the two-leg ladders introduced above 
for different numbers of sites and different values of the coupling strengths 
Jt, Ji and Jc in the Mtot = subspace. 

First test: L= 9, Jt=15 and 2.5, Ji=5, Jc=3 

Results for L = 9 sites along a leg in the Hilbert space spanned by 
= 48620 basis states with Mtot = are presented in Figs.(l2](a-f)): 
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• For Jt = 15, Fig.([2](a)) shows that p{l) increases from very small values 
to 0.5% when n lies between 500 and 80, p{2) and p{3) ~ 0.05% when 
n ~ 500. 



• For Jf = 2.5, Fig. ([21^c)) shows that p{l) increases also from very small 
values to 0.5% when n lies between 1000 and 100, p(2) and p{3) ~ 2.5% 
when n ~ 500. 

Instabilities in the position of the eigenstates are localized at dimensions 
of the reduced Hilbert space where the elimination of components of the 
wave function cannot be properly compensated by the renormalization of 
g = Jt which changes sizably at these values. This can be shown by the close 
relation in reduced Hilbert space between the behaviour of the renormalized 
coupling parameter J^"^ and the entropy s which works as a measure of the 
distribution of the amplitudes in the ground state, see Figs.([2]^b)-12]|^e)) and 
(l2](d)-[2](f)). The effect of instability is weaker for Jt = 15 which corresponds 
to a strong coupling along the rungs whereas Jt = 2.5 corresponds to a 
stronger coupling along the legs of the ladder. This can be explained by 
means of symmetry arguments [HE]. 

Second test: L= 6, Jt=15, Ji ^ 12.21, J^ ~ 12.21, 11 

Here we test the behaviour of the system at a first order phase transition 
point. 

For L = 6 sites along a leg the Hilbert space is spanned by = 924 basis 
states with Mtot = 0. Results are shown in Figs([3](a-d)). 

This case corresponds to the transition from a rung dimer phase to a 
Haldane phase which appears for J/ = Jc when {Jt/ Ji)crit — 1.401 in the case 
of an asymptotically large system [TS]. The ratio depends on the size of the 
system. As predicted by the theory above the coupling constant g = Jt is 
expected to stay constant at the level crossing point. 

For Jt=l5, Ji = Jc^ 12.21, ci = -11.25,62 = -11.25,63 - -10.7. One 
observes a level crossing between the ground state and first excited state of 
the energy spectrum. Fig.([3](b)) shows the constancy of the coupling strength 
down to low dimensions of < 100. For smaller N, when components of 
the wave function with sizable weight are eliminated, Jt increases sizably. 
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In Fig.([2]^a)) one notices that the ground state remains stable but the first 
excited states move abruptly during the reduction procedure. 

For Jt = 15, Ji - 12.21, J, ~ 11, d - -11.38,62 ~ -10.67,63 - -10.4 
one stays in the vicinity of the transition point. Figs.([21^c-d)) show the be- 
haviour of the spectrum and the coupling strength Jt. One notices that in 
Fig.([3t^d)) the renormalized coupling parameter stays stable down to < 
250. In Fig.([3](c)) the ground state remains stable but the first excited states 
goes moving during the reduction procedure but less than in the case shown 
in Fig.p:a)). 

Summary - conclusions. 

In the present work we tested and analysed the outcome of an algorithm 
which aims to reduce the dimensions of the Hilbert space of states describing 
strongly interacting systems. The reduction induces the renormalization of 
the coupling strengths which enter the Hamiltonians. By construction the 
algorithm works in any space dimension and may be applied to the study of 
any microscopic A^-body quantum system. 

The outcome of numerical tests of the method applied to strongly corre- 
lated and frustrated quantum spin ladders can be summarized as follows: 

• Local spectral instabilities of the ground and low excited states ap- 
pearing in the course of the reduction procedure developed above are 
correlated with the elimination of basis states with sizable amplitudes 
in the ground state wavefunction. The renormalization is able to cure 
this instability down to a small values of the Hilbert space dimensions. 
Numerical examples show that the procedure allows for sizable reduc- 
tion of the dimensions of Hilbert space. 

• The stability of the low-lying states of the spectrum in the course of 
the reduction procedure depends on the relative strengths of the cou- 
pling constants. Here the ladder favours a dimer structure (i.e. strong 
coupling along the rungs) for which the stability is the better the larger 
this coupling. Symmetry arguments can explain such a behaviour. 

• As it could be expected the system is unstable at phase transition 
points. The theory predicts a constant coupling parameter there. This 
is the case in numerical applications, down to a limit related to the 
space dimension reduction discussed above. The ground state stays 
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remarkably stable, the spectrum of excited states gets however strongly 
unstable there. 

Last, the present procedure may be extended to the renormalization of 
more than one coupling parameter and possibly to systems which are not 
linear in g. However in most cases realistic microscopic systems and in par- 
ticular quantum spin systems are described by Hamiltonians which depend 
linearly on coupling strengths. The study of systems at finite temperature 
can also be performed ^24j. 
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Figure 2: Behaviour of {pi,i = 1,2,3}, Jt and the entropy s as a function of 
the size n of the Hilbert space. Cases (a) — (6) and (e) correspond to Jt = 15, 
(c) — (d) and (/) to Jt = 2.5, J; = 5, Jc = 3. The number of sites along a leg 
is L = 9. 
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Figure 3: Same notation as in Fig(j 



Jt 



Jt 



20 
19 
18 
17 
16 

15 L 



+ 
+ 








1 


Jt 


^ 






1 


















+ 


















+ 
+ 

+ 
+ 



























































100 200 300 400 500 600 700 800 900 
n (b) 



20 
19 
18 
17 
16 
15 



T r 



+ 
















+ 
































i 

















T 1 1 r 

Jt + 



(Jt/Jl) 



crit 



1.23, Jt 



15, Ji 



Jr 



100 200 300 400 500 600 700 800 900 
n (d) 

|. Cases (a) — (6) correspond to 
12.21. Cases (c) — (d) correspond 



to Jt = 15, Ji ~ 12.21, Jc = 11. The number of sites along a leg is L 
Broadened lines are drawn in order to facilitate the reading. 



6. 



13 



